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Abstract 

The well-known "displace, cut, and reflect" method used to generate cold 
disks from given solutions of Einstein equations is extended to solutions of 
Einstein-Maxwell equations. Four exact solutions of the these last equations 
are used to construct models of hot disks with surface density, azimuthal 
pressure, and azimuthal current. The solutions are closely related to Kerr, 
Taub-NUT, Lynden-Bell-Pinault and to a one-soliton solution. We find that 
the presence of the magnetic field can change in a nontrivial way the different 
properties of the disks. In particular, the pure general relativistic instability 
studied by Bicak, Lynden-Bell and Katz [Phys. Rev. D47, 4334, 1993] 
can be enhanced or cured by different distributions of currents inside the 
disk. These currents, outside the disk, generate a variety of axial symmetric 
magnetic fields. As far as we know these are the first models of hot disks 
studied in the context of general relativity. 
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1 Introduction 



The existent magnetic fields in the primeval cloud of plasma that originates 
stellar objects like galaxies and stars remain after the formation of the ob- 
jects. The magnetic lines of force are said to be "frozen in" the fluid. This 
field can play many different roles in astrophysics. Two examples are: The 
radius of a white dwarf for extreme relativistic degeneracy can change signif- 
icantly even when a small quantity of magnetic energy is present . Binary 
and active-galactic-nuclei disks are hot enough to be significantly or fully 
ionized, and their high conductivity allows the current to flow freely. The 
presence of magnetic fields strongly influence the structure and evolution of 
these disks p. 

Axially symmetric solutions of Einstein field equations corresponding to 
disk like configurations of matter are of great astrophysical interest and have 
been extensively studied. These solutions can be static or stationary and 
with or without radial pressure. Solutions for static disks without radial 
pressure were first studied by Bonnor and Sackfield and Morgan and 
Morgan and with radial pressure by Morgan and Morgan ||, and, in 
connection with gravitational collapse, by Chamorro, Gregory and Stewart 
IJ. Disks with radial tension has been recently considered in 0. 

The stability of models with no radial pressure can be explained by 
either assuming the existence of hoop stresses or that the particles on the 
disk plane move under the action of their own gravitational field in such 
a way that as many particles move clockwise as counterclockwise. This 
last interpretation is frequently made since it can be invoked to mimic true 
rotational effects. A large class of thin disks solutions were obtained by 
Letelier and Oliveira ||] . Solutions for self-similar static disks were analyzed 
by Lynden-Bell and Pineault ||, and Lemos |10|j . The superposition of static 
disks with black holes were considered by Lemos and Letelier [11, 12, [HJ, 
and Klein [14|. Recently Bicak, Lynden-Bell and Katz [15| studied static 
disks as sources of known vacuum spacetimes and Bicak, Lynden-Bell and 
Pichon ]l6| found an infinity number of new static solutions. 

The above mentioned solutions of the Einstein equations have no sources 
other than the disk itself, i.e., outside the disk their metrics are solutions of 
the vacuum Einstein equations for static spacetimes with axial symmetry. 
The search of new solutions is facilitated in this case by the fact that one of 
the Einstein equations is equivalent to the usual (linear) Laplace equation. 

In the present paper we study models of disks that are characterized by 
a surface density, an azimuthal pressure and a current density. Out side the 
disk the metric satisfy the Einstein-Maxwell equations for a static axially 
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symmetric spacetime filled with a magnetic field. As far as we know these 
are the first models of hot disks studied in the context of general relativity. 

We find that the presence of the magnetic field can change in a nontriv- 
ial way the different properties of the disks. In particular, the pure general 
relativistic instability studied by Bicak, Lynden-Bell and Katz [JD]] can be 
enhanced of cured by different distributions of currents that outside the disk 
generate magnetic fields of different forms. 

In Sec. II we study the Einstein-Maxwell equations for a static axially 
symmetric spacetime. We introduce the well known "displace, cut, and 
reflect" method to generate a disk from a given solution of these equations. 
We also present the expressions for the main physical variables of the disk. 
In the next section, Sec. Ill, four models of disks are studied. They are 
based in simple solutions to the Einstein-Maxwell equations. The first two 
are similar to Kerr and Taub-NUT solutions. The third is a generalization 
of a Lynden-Bell and Pinault solution Q and the last is one obtained using 



solitonic techniques [17]. In the last section , Sec. IV, we discuss some of 



the found results and make some comments. 



2 Einstein-Maxwell Equations and Disks 

The simplest metric to describe static and axially symmetric systems is the 
Weyl metric, 

ds 2 = fdt 2 - r l [e w {dz 2 + dr 2 ) + r 2 d<p 2 ], (1) 

where the functions / and w depend only on the cylindrical coordinates r 
and z. 

The Einstein-Maxwell system of equations is equivalent with, 



G p.v = _ T ^ ; ( 2 ) 

T „u = F P aF a» _ lg^F a PF a/3 , (3) 

V a F a/1 = 0, (4) 

F^u = d»A v - cM M , (5) 

where all the symbols have their usual meaning. Now we choose the magnetic 
potential as 

A^A^zW, (6) 
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in other words we take the electromagnetic potential for a pure axially sym- 
metric magnetic field, 



F ipg = -d s A, F rv = d r A. (7) 

Therefore the physical components of the magnetic field are, 

Br = ~(f/r)e~ w/2 d z A, (8) 
B z = {f/r)e~ w l 2 d T A. (9) 

From (p])-(||) we find the well known system of Einstein-Maxwell equations 
for axially symmetric spacetimes with pure magnetic fields [18|, 



r 2/2 + r ' ( 0) 

, rd r fd z f 2fd r A8 z A 

d z w = /2 + , (11) 

d r (rd r f/f - fAd r A/r) + 8 z (rd z f/f - fA8 z A/r) = 0, (12) 
d r (fd r A/r) + 8 z (fd z A/r) = 0. (13) 

Once the solution of the system of equations (|l2|)-(|l3|) is known, Eqs. ( |l0| ) 
(|IT|) give us the function w as a quadrature whose existence is guaranteed 
by the previous equations. 

The method that we shall use to generate the metric of the disk, as 
well as, its material and electromagnetic content will be the well known 
"displace, cut and reflect" method that was first used in Newtonian gravity 
by Kusmin and Toomre [^(J and later in general relativity by most of 
the authors that write about disks and also fine shells, defects, etc. 

The material content of the disk will be described by functions that 
are distributions with support on the disk. The method that we shall use 
can be divided in the following steps: First, choose a surface that divide the 
usual space in two parts; one with no singularities or sources and the other 
with the sources. Second, disregard the part of the space with singularities. 
Third, use the surface to make and inversion of the nonsingular part of the 
space. The result will be a space with a singularity that is a delta function 
with support on the surface. This procedure is depicted in Fig. 1 for the 
gravitational field of a punctual mass. 

In the present case the procedure is equivalent to make the transforma- 
tion z — > zq + \z\, with zq constant. In the Einstein equation for axial sym- 
metry we have first and second derivatives of z. Note that d z \z\ = 29{z) — 1 
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and d z 0(z) = 5(z) where 9{z) is the Heaveside function and 5{z) the usual 
Dirac distribution. 

The procedure above mentioned applied to the metric (|l]) and then 
Einstein-Maxwell equations gives us 

G^ = -(T^ + f^), (14) 

V M F^ = r, (15) 

with 

f° = p = e- w (2d z f-fd z w)5(z), (16) 

f$ = P< p = e- w fd z w6(z), (17) 

J* = 2fd z A5(z)/r, (18) 

where d z f should be understood as d z f\ z= o+, etc. Tf*™ is the electromag- 
netic energy-momentum tensor defined in (|3|), T^ u is the matter energy- 
momentum tensor on the plane z = and J M is the current density on the 
plane z = 0. The "physical meassure" of length in the direction d z for the 
metric ([]]) is (e w / fY^dz, then the invariant distribution is 5(z)/(e w / /) 1//2 . 
Therefore the surface density and the planar pressure or tension along the 
direction d v are, 

° = [{e' w /f) l /\2d z f - fd z w)] z=0+ , (19) 
P^ = [(fe- w ) 1/2 d z w] z=0+ . (20) 

In the same form one can compute the tetradical component of the planar 
current density, we find, 



2[e w ' 2 d z A] z=Q+ . (21) 



Another physically meaningful quantity associated to non rotating disks is 
the counter-rotation velocity, 



r2 l_ f d zW 



2d z f - fd z w 



(22) 

=0+ 



The disk is supposed to be formed by the same number of particles moving in 
circles in one direction than in the opposite direction. The specific angular- 
momentum of a particle with rest mass \i rotating at a radius r, defined as 
h = -Kyi ix = g^dtp/ds, is 



h = T f^ . (23) 



a 2 -p* 



5 



This last quantity can be used to detect instability of the disk against ra- 
dial perturbations (ring formation). The condition of stability is usually 
presented as 

V > °- (24) 

i.e, for h > we have stability when the specific angular momentum is an 
increasing function of r for < r < R (R the radius of the disk). This 
criteria is an extension of Rayleigh criteria of stability of a fluid at rest in a 
gravitational field, see for instance [21] . For other significant quantities for 
relativistic disks sec [15]. 

As we said before in the present case the spacetime is static, hence in 
the disk we do not have a net movement of matter in one direction, but 
we have current along the cL, direction. This means that we have a kind of 
moving massless charge approximation. 

The field lines of the magnetic field are given by the ordinary differential 
equation, 

t = t- (25) 

From (||) and @ we find, 

d r Adr + d z Adz = 0. (26) 

Thus the equation A(r, z) = C (C constant) gives as the lines of force of the 
magnetic field. 



3 Exact Disks with Magnetic Fields 

There are many techniques to find exact solutions of fll2"|)— fll3|), they are 
based in the fact that this system of equations has extra symmetries (Geroch 
group p2|) and that is closely related to the principal sigma model |L7]]. In 
particular, we can write (|l2])-([l3|) in the matrix form 

d r [rd r (M)M- 1 )} + d z [rd z (M)M~ 1 )} = 0, (27) 

where 

M= [ (f/2)WA 2r2 + fA2 )' (28) 

d r [™ + ln(r 4 // 2 )] = -Tr[((d r M)M- 1 ) 2 - ((d^M- 1 ) 2 }, (29) 

r 

OJw + ln(r 4 // 2 )] = -Tr\(d r M)M- 1 )(d z M)M- 1 )}. (30) 

r 



2(/) 



1/2 



as well as (0)-(0), 
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Note that the matrix M is symmetric and det(M) = r 2 . Eq. (27) is mathe- 
matical equivalent to the Ernst equation, although the physical and geomet- 
ric content of these equations are different Jlq] . Now we shall study models 
of disks based in simple solutions of (p?|) 

A. Disk Associated with a Magnetic Dipole Solution 

A Kerr like solution to (|27]) is 

,1/2 p 2 x 2 - g 2 y 2 - 1 r „ n 

1 (px+l) 2 -<7V' 1 J 

A = ^p /^rf + 2 f 7 2) , (32) 
p z ar — q z y z — 1 

where x and y are the spheroidal coordinates 



2px = Jr 2 + (z+p) 2 + Jr 2 + (z-p) 2 , (33) 



2py = sjr 2 + (z+p) 2 - ^Jr 2 + (z-p) 2 , (34) 
p and g are two constants related by 

q 2 = p 2 - i- (35) 



The function w is given by 
with 



p 2 x 2 — q 2 y 2 — 1 



p 2 ((px + l) 2 — q 2 y 2 ) ^ 
A similar solution, but in a different system of coordinates, was first studied 



by Bonnor [23], for p = 1 reduces to the Weyl 7-metric (with 7 = 2). 
This solution can be interpreted as a Weyl 7-metric with a magnetic dipole 
located in the origin of the coordinate system. 

Now we perform the "displace, cut and reflect" method above men- 
tioned, i.e., we do z — > zq + \z\ in the Kerr like solution. This is equivalent 
to change only the definitions of the coordinates x and y by the new defini- 
tion 



2px = ^r 2 + (|*| +z +p) 2 + yr 2 + (|z| +z -p) 2 , (38) 
2py = \jr 2 + (\z\ +z +p) 2 - ^r 2 + (\z\ + z -p) 2 . (39) 
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Now, the computation of the physical quantities like cr,p^, etc. associ- 
ated to this disk is straight forward, but the result is rather cumbersome, 

a = 8[({2yz - 3)x 2 - y 2 - (x 2 + y 2 )px)p 2 y + (x 2 + y 2 - 2q 2 x 2 y 2 )z 
+ {(3yz - 2)y + x 2 z + (x 2 + f)q 2 y)px]p 2 / 

\{q 2 y 2 + 1 -p 2 x 2 )(qy + 1 +px) 2 {qy - 1 - px)\ (40) 
p v = 8(((x 2 + y 2 )q 2 + 2)pxy - (x 2 + f + 2q 2 x 2 y 2 )z - ({x 2 + f)p - 
2xyz )p 2 xy)p 2 /[(q 2 y 2 + 1 - p 2 x 2 ) 2 (qy + 1 +px)(qy - 1 - px)], 

(41) 

j v = -4^[(((3y 2 - l)x 2 z + 4y + 2(x 2 + f)q 2 y)px - 
({y 2 + 2yz - l)y - 2(y - z )x 2 + (4x 2 yz - 2x 2 + 

y A - y 2 )q 2 y))p - {(y 2 + 1 - (y 2 - 3)? 2 y 2 )zo + 2(x 2 + y 2 )/y 

+(y 2 - Ayz - 1 + 2x 2 )p 3 xy)x]g/(x 2 - y 2 ) 3 , (42) 

where 

2px = yV 2 + (z +p) 2 + ^/r 2 + (z -p) 2 , (43) 

2py = ^r 2 + (z +p) 2 - \Jr 2 + (z -p) 2 . (44) 

The magnetic potential A is given by the same function (|32] ) but with x and 
y given by © and ©. 

We shall study the above functions graphically. In Fig. 2 we present 
the surface density a as a function of the radial coordinate r for zq = 2 and 
different values of the parameter q that controls the magnetic field. 

The density have a finite central maximum that decays fast. The curve 
with the lowest maximum has q = and the next three curves have q = 
0.5, 1, and 1.5. The curve with q = represents the density of a disk with 
no magnetic field, i.e., the density of a usual cold disk. We see that the 
inclusion of the magnetic field increases the central density for the values of 
q shown. Really, the central density increases up to a value close to 2.5 for 
q = 1.73 and for values of q > 1.73 decreases very fast and take negative 
values (not shown in the figure). 

The azimuthal pressure is shown in Fig. 3 for zq = 2 and q = (bottom 
curve), 0.5,1,1.5 (top curve). We have that the pressure has a maximum 
that increases with q and moves toward the disk centre. For g > 1.5 we have 
similar behavior. 

The azimuthal current is presented in Fig. 4 for the same values of the 
parameters zq and q . The curve that corresponds to q = is the axis r. The 
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current has similar behavior than the pressure. We see that the behavior 
of the current and the pressure is consistent with the presence of a ring-like 
current of massless particles. 

In Fig. 5 we present the magnetic potential A and its level curves that 
represent the the magnetic field for zq = 2 and q = 1. We see that the form 
of the magnetic field also supports the idea of a ring-like current distribution. 

Finally for this disk, we shall examine the stability against radial per- 
turbations (ring formation). In Fig. 6 we present the radial dependence 
of the angular momentum h for a cold disk q = for zq = 2 (top curve), 
zq = 2.5,4,3.4 and 4 (bottom curve). We see the the disk with zq = 2 is 
unstable, h is not an increasing function of r for < r < 14. This is a pure 
general relativistic instability. Now we shall study the stability of the disk 
with zq = 2.5 when we switch on the magnetic field. 

The radial dependence of the angular momentum h is presented in Fig. 
7 for q = (bottom curve), 0.5, 1, 1.5 (top curve). Therefore the disk with 
q = 1.5 is unstable. Again the magnetic field introduces a pure gravitational 
instability. Note that the disk with q = is stable and the inclusion of the 
magnetic field makes it gravitationally unstable. The density, pressure and 
current for this disk look similar than the ones for zq = 2, we shall comeback 
to this point later. 



B. A disk associated with Taub-NUT like charge 



A Taub-NUT like solution is, 

/2 = pV-i) ( 

px z + Ix + p 
A = 2V2by/p, (46) 

F = PX ] + 2 2X+ 2 " (47) 

with 

b 2 = l-p 2 . (48) 

The function w is given by (|36"|), b is a constant and x and y are the same 
coordinates defined before. This solution for p = 1 also reduces to the Weyl 
7- metric (with 7 = 2). 

As in the precedent case we find, 

° = "8[((y 2 + 2 + x 2 )p 2 y - (x 2 + 3y 2 )z )x - ({2y 2 z - 3y + 
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z )x 2 - (y - z )y 2 )p]/[((x 2 + l)p + 2x) 2 (x 2 - l)p], (49) 
-8[((2y 2 - l)x 2 - y 2 )z -(f-2 + x 2 )pxy] 
Plp ((x 2 + l)p + 2x)(x 2 -Ifp ' 1 J 

j v = AV2(px-yz )(x 2 -l) 2 b/[(x 2 -f) 3 p 2 ]. (51) 

The magnetic potential A is the given by the same function (^) but with 
x and y given by (pq) and (|39|). 

The surface density as a function of r is presented in Fig. 8 for 6 = 
(bottom curve), 6 = 0.3, 0.6, 0.9 (top curve) and zq = 2. We see that the the 
density maximum increases when the magnetic parameter b increases. 

We also see (Fig. 9), for the same value of the parameters, a similar 
behavior for the current. The curve b = is the axis r. The pressure behaves 
in the opposite way (Fig. 10), it decreases when the magnetic parameter 
b increases. For the top and bottom curves we have 6 = and 6 = 0.9 
respectively. 

The field lines of magnetic field in this case are quite simple, they are 
just parts of hyperbola branches that cut the disk in an axially symmetric 
way, the axis of the disk being the symmetry axis. Note that the distribution 
of current is quite different in this case than in the precedent one. Now we 
have a more homogeneous distribution. 

The angular momentum h as a function of r is shown in Fig. 11 for 
zq=2 and 6 = (top curve), 6 = 0.3,0.6,0.9 (bottom curve). We see that 
the inclusion of the magnetic parameter can stabilize the disk. We shall 
comeback to this point later. 

C. A generalization of a Lynden-Bell-Pinault solution 



A magnetized Lynden Bell-Pinault like solution is 

/ 1/2 = a 2 //o + cV//o, (52) 

A = -V2 ah {° + d 2 i (53) 

w = «j + 21n(///o)-81na, (54) 

where /o and i^o are the respective functions associated to the Lynden-Bell- 
Pinault solution 101: 



f = [( z+ V^T^)/2] 2n , (55) 

2 z + Vr 2 + z 2 
w = An In— —==-, (56) 
2v r + z 
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a, b, c, d and n are constants restricted by the relations ad—cb = ±1, and n < 
1. This solution was found by applying the transformation M = B _1 MqB 
with 

B =(:i)- <") 

Mo being the matrix fl2q ) com A = 0. Note the right hand side of Eqs. ( |29| ) 
and ( p0| ) are invariant under the transformation M = B~ 1 MqB. 

Now doing z — > |z| + i n the previous solution we find after doing 
z = 0, 



a = r 3« 2 2n 2 -n +2(1 _ ^4^/ [(r 2n fl 2 + 2 2n c 2 r 2 )r 



r 2n a 2 n — 2 2n (2 — n)c 2 r 2 ]a A n 



j s _2 2n + 5 l 2 - 



■2n' 



(j-2n a 2 _|_ 22n c 2 r 2^2 
l-2n /„3 



cnr 



(58) 
(59) 
(60) 



Note that this solution has well behave density that is singular in the center 
of the disk, like the original LBP solution. The pressure is positive in the 
center, but for r > r c it becomes negative, 



na 



2 2n (2 -n)c 2 



l/(2-2n) 



(61) 



i.e., we have tensions. This is a general feature of solutions generated in the 
same form using a Weyl solution as a seed. 

Hence in this case we have that this particular magnetization of the 
LBP disk introduces tensions. 



D. A particular disk associated to a 1-soliton like solution 
A particular 1-soliton solution is 



/ 



1/2 



(eg + c\)iir 



C 2 fl 2 



c\r 2 ' 



A 



F 



V2- 



c qCiQu 2 +r 2 ) 

(co + c?> 
CgU 2 - c\r 2 

\/r(c 2 + c 2 )(n 2 + r 2 ) ' 



(62) 
(63) 
(64) 
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where 

\x = z - zq + \jr 2 + {z- z ) 2 (65) 

and cojCi and zo are constants. This solution was obtained by applying 
the Belinsky-Zhakarov [o] method with a single pole scattering matrix. A 



similar solution was studied in [24] 



As before by doing z — > \z\ we find 



a = 16(cg + c?)/i 2 r 2 /[(c ^ + cir) 2 (coM-cir) 2 ], (66) 
£V = 8 (^ 2 - r2 ))/[( c o^ + cir)(c /i-cir)] (67) 
j v = AV2rix 2 {i? -r 2 ) Co ci/(/i 2 + r 2 ) 3 (68) 

For this solution we have that the central density of the disk is zero but not 
the pressure. 

In the four solutions presented we have chosen the integrations con- 
stants associated with the quadrature of w in such a way that w — > on the 
symmetry axis. Thus we have a regular symmetry axis out side the sources. 
There are many known exact solutions of the Einstein-Maxwell equations 
with axial symmetry, as well as, many methods to find them. In particular, 
the first two solutions presented are specializations of a Kerr-NUT like so- 
lution or two soliton solution. We have presented them separated because 
of algebraic simplicity. 



4 Discussion 



For the disks studied in the present paper we made a rather complete study 
of parameters, we only mentioned the more relevant aspects of each solution. 
In particular, for solutions generated from finite sources like the one in 
Fig. 1, we have that the density of the disk is a decreasing function of the 
"cut" parameter zq. For a Newtonian disk generated from Newton's law 
$ ~ M/\/r 2 + z 2 the "displace, cut, and reflect method" gives the surface 
density ~ Mzo/(r 2 + Zq) 3 / 2 (Kusmin-Toomre disk). Thus the central 
density o"at(0) ~ M/zq. For the first two studied disks A. and B. we put 
the " mass" parameter equal to one. This is equivalent to have coordinates 
t, r, z dimensionless, one can recuperate this parameter by multiplying the 
line element by m. The numerical study of the two first disks shows that 
the central density also decays with zq, but with a different law than the 
Kusmin-Toomre disk. 
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We find that the magnetic field stabilize the disk against radial pertur- 
bation when the distribution of current that generate the field has its origin 
in a NUT like parameter. And made the disk unstable when the distribu- 
tion of current is originated by a Kerr like parameter. In a different, context 
we studied a similar phenomena, we found that the inclusion of magnetic 
mass in a Weyl solution with unstable orbits (Schwarzschild solution plus 
a dipolar halo) makes the orbits less chaotic p5| , i.e., more stable. Also 
the opposite is true for particles moving in counter rotation in a slow rotat- 
ing Kerr solution plus a dipolar halo. The rotation makes the orbits more 
chaotic pq| . These are pure general relativistic effects. 

In the present work we studied simple four solutions of the Einstein- 
Maxwell equation generated by well known new-solution-generating algo- 



rithms [18|. There are quite a few solutions of the Einstein Maxwell equa- 
tions with axial symmetry that can be used to generate disks. For instance 
in the present context members of a Tomimatsu-Sato like family of solutions 
are good candidates. But the algebraic complexity involved in the finding 
of the physical variables of the disks is not trivial. 

The hot disks studied in this paper have not radial pressure nor the 
charges have mass. There are several methods to add radial pressure to 
disks, they are based in the fact that a coordinate transformation that keep 
detM as a harmonic function in the variables r, z does not change the form 
of the metric (|l|). 

A second degree of complexity is to add rotation to the disks, also suit- 
able solutions of the Einstein-Maxwell equations for stationary spacetimes 
are well known. Again the algebraic complexity is not trivial in this case. 
The inclusion of rotations, in principle, can be used to introduce mass in 
the charge carriers. A particular case of hot rotating disk is under study. 

Finally, we want to mention that all the computation of this work was 



performed using the algebraic programming system Reduce |27] 
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Fig. 1. The illustration of the "displace, cut, and reflect" method for the 
generation of a disk from the gravitational field of a punctual mass. In A we 
displace and cut out the singularity with a plane, in B we reflect the field 
on the plan. We end up with the field of a disk. 
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Fig. 2. The surface density a as a function of r for zo = 2 for q = 
(lower maximum), 0.5, 1, 1.5 (higher maximum). The parameter q controls 
the magnetic field. 
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Fig. 3. The azimuthal pressure for zq = 2 and q = (bottom curve), 
0.5,1,1.5 (top curve). We have that the pressure has a maximum that 
increases with q and moves toward the disk centre. 
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Fig. 4. The azimuthal current for the same values of the parameters z$ and 
q. The curve that corresponds to q = is the axis r. 
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Fig. 5. The function A and its level curves that represent the the magnetic 
field for zq = 2 and q = 1. We see that the form of the magnetic field 
support the idea of a ring-like current distribution. 
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Fig. 6. The radial dependence of the angular momentum h for a cold disk 
q = for zq = 2 (top curve), zq = 2.5, 4, 3.4 and 4 (bottom curve). 
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Fig. 7. The radial dependence of the angular momentum h for q = 
(bottom curve), 0.5,1,1.5 (top curve). In this case colder disks are more 
stable. 
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Fig. 8. The surface density as a function of r for 6 = (bottom curve), 
b = 0.3, 0.6, 0.9 (top curve) and zq = 2. The density maximum increases 
when the magnetic parameter b increases. 
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Fig. 9. The current for the same value of the parameters of Fig 8. The curve 
b = is the axis r. The current and the density have a similar behavior. 
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Fig 10. The pressure decreases when the magnetic parameter b increases. 
For the top and bottom curves we have 6 = and b = 0.9, respectively. 
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Fig. 11. h as a function of r for zq=2 and b = 0(top curve), b = 0.3, 0.6, 0.9(bottom 
curve). In this case hoter disks are more stable. 
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